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1. Introduction 



Our modified gravity (MOG) theory, also known in previous work [T] as Scalar- Tensor- 
Vector Gravity or STVG, is based on an action that incorporates, in addition to the 
Einstcin-Hilbert term and the matter action, a massive vector field, and three scalar 
fields corresponding to running values of the gravitational constant, the vector field 
coupling constant, and the vector field mass. 

A Dirac-Hamiltonian analysis [5] can be used to show that the theory is free 
of constraints and it is free of ghosts and instabilities. The theory has been used 
successfully to explain cosmological observations [3J, the motion of galaxy clusters 
[I] , the Bullet Cluster [5] , the rotation curves of galaxies and dwarf galaxies [5] , and 
the velocity dispersion profiles of satellite galaxies [3 H] and globular clusters [9] 
without exotic dark matter. However, until now these explanations required case-by- 
case fitting of two parameters, or alternatively, ad-hoc formulae, not derivable from 
the action principle, to predict the values of these parameters. 

In this paper, we demonstrate that it is in fact possible to eliminate fitted 
parameters altogether in favor of universal constants. After fixing initial conditions 
using observational data, the field equations can be integrated numerically, and 
the theory can be used to make predictions from the scale of the solar system to 
cosmological scales, i.e., across at least 14 orders of magnitude in length, or more than 
20 orders of magnitude in mass-energy. 

We begin in Section [2] by introducing the theory through the action principle, 
and establish key assumptions that allow us to analyze physically relevant scenarios. 
In Section[3j we derive the field equations using the variational principle. In Section|4] 
we solve the field equations in the static, spherically symmetric case. In Section [5l we 
postulate the action for a test particle, and obtain approximate solutions to the field 
equations for a spherically symmetric gravitational field. In Section [5] we demonstrate 
how the Friedmann equations of cosmology can be obtained from the theory. In 
Section [3 we utilize the theory to obtain new estimates for galaxy rotation curves, 
derive the Tully-Fisher law, and show how the solutions we obtained for the field 
equations remain valid from cosmological to solar system scales. Lastly, we end in 
Section [5] with conclusions. 



2. Modified Gravity Theory 



The action of our theory is constructed as follows pQ. We start with the Einstein- 
Hilbert Lagrangian density that describes the geometry of spacetime: 

C G = - T J- 7 ^(R+2A)^, (1) 



16ttG 

where G is the gravitational constant, g is the determinant of the metric tensor 
(we are using the metric signature (+, — , — , — )), and A is the cosmological constant. 
We set the speed of light, c = 1. The Ricci-tensor is defined as 



Ruv 



"a 1 [iv 



U V L HO 
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where T" is the Christoffel-symbol, while R ■ 
We introduce a "fifth force" vector field 
density: 



pa -p/3 



(2) 



g J V' 

bp via the Maxwell-Proca Lagrangian 



v Air 



(3) 
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where £? M „ = d^v—d^^, fi is the mass of the vector field, u> characterizes the strength 
of the coupling between the "fifth force" and matter, and is a self-interaction 
potential. 

Next, we promote the three constants of the theory, G, [i and u>, to scalar fields 
by introducing associated kinetic and potential terms in the Lagrangian density: 



C s = - 



G 



G 2 



^9, 



(4) 



where V M denotes covariant differentiation with respect to the metric g llUl while Vq, 
V^, and Vu are the self- interaction potentials associated with the scalar fields. 
Our action integral takes the form 



(C G + £<{, + C s + C M ) d 4 x, 



(5) 



where Cm is the ordinary matter Lagrangian density, such that the energy-momentum 
tensor of matter takes the form: 

2 SSm 



T — — - 



-9 Sg^ 

where Sm = J Cm d 4 x. A "fifth force" matter current can be defined as: 
j V _ 1 SSm 



(6) 



(7) 



We assume that the variation of the matter action with respect to the scalar fields 
vanishes: 

SS M 



SX 

where X = G,fi, u). 



= 0, 



(8) 



3. Field equations 



The field equations of the theory can be obtained in the form of the first and second- 
order Euler-Lagrange equations: 

w 



— - V 

dX v 'd(9 v X) 



dC 
dC 



dC 



0(V„« 



= o, 



d K + d K d\- 



dC 



L fivi 



(10) 
(11) 



dg^ K dg^, K "dgi»', K x 

where C — Cq + C^ + Cs + Cm is the total Lagrangian density and a comma in the 
covariant index is used to indicate a coordinate derivative. 
The full set of the theory's field equations reads: 



1 

47T 



V>V„m V v GV„n 1 



(12) 



// 



G 



47T 



G^^r+-vM-y'M = o, (13) 
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+K» = o, 



(14) 



V'V.G- 



-V^(G) + G 
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G 



A* 



+ — (i* + 2A) = 0, 
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47T 
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+ ^(^gv^) + ^1 + YM + vUu 

+R^ - ^g^R + g^A = -8nGT„„. 



(15) 



(16) 



4. Static, Spherically Symmetric Vacuum Solution 



In the static, spherically symmetric case with line element 

ds 2 = Bdt 2 - Adr 2 -rW, 
with dQ 2 = d9 2 + sin 2 0d(f> 2 , the field equations are written as 



dV 6 



47T 



(17) 
(18) 



J t + - ft 
r 
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B 



b' t - Afi 2 4> t + AB 



d<p t 



An A 



CO 



(19) 



+AV^{G) - 3A^§1 - \G 



G 



m 2 



2 V B A. 

AG(R + 2A) 



+ V u (u) 



16tt 



= 0, 



(20) 
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(21) 
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where 



J ' + r J -^' + ^\B^~^) + 2KB A 
+ 8^B G ^ -L AGV ^ AV ^ = °' (22) 
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,u 2 G 2 
V = 2uGV 4> (<p) + 8tt 



G 2 fi 2 



(29) 



The prime denotes differentiation with respect to r, i.e., y' = dy/dr. 

These equations can be substantially simplified in the matter vacuum case 
(Tjf = 0), with no cosmological constant (A = 0), setting the potentials to zero 
(V<f, = Vg = Vfx = Vui = 0) and also setting <j> r = 0. These choices leave us with six 
equations in the six unknown functions A, B, <f> t , G, fj,, and lj, which read after some 
trivial rearranging: 



2ABG AGr \AB B ) 2AB 4AB 2 

A'B' B' A' 1 1 
+ 4A 2 B ~ 2AB~r ~ 2AV + Ar 2 ~ r 2 ~ °' (30) 

G" 2G' 2 B'G' A'G' B' A' 

AG ~ AG 2 ~ 2ABG ~ 2A 2 G + ~ABr + AV 

, / G' 2 ^' 2 J 2 LuGfi 2 
Btt — — + -^7 ; A 
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Figure 1. Comparing MOG numerical solutions to the Reissner-Nordstrom 
solution, for a 10 11 Mq source mass. The MOG parameters A (solid red line) 
and B (dashed brown line) are plotted along with the Reissner-Nordstrom values 
of A (dash-dot green line) and B (dotted blue line), calculated using 11511 . 
Horizontal axis is in pc. We observe that the A parameter reaches at below 
the Schwarzschild radius of a 10 11 Mq mass, which is ~ 0.01 pc. 
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(32) 
(33) 

0, (34) 
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r ~2~G 
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(35) 



A solution to this set of six equations can be obtained using numerical methods. 
This solution is exact in the sense that it is obtained without dropping any terms, 
and its accuracy is limited only by the machine precision and the stability of the 
integration algorithm. This approach requires that we find suitable initial conditions 
for the six unknown functions A, B, G, fi, to, and <p t . This is accomplished in the next 
section, using the motion of a point test particle as a guide. 

Before we proceed to the case of the test particle, however, let us describe in 
approximate terms the type of solution that we obtain through numerical methods. 
The following findings apply when the distance r from the source mass is much larger 
than the Schwarzschild radius r$ = 2GM associated with the source mass, r r$- 
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All three scalar fields turn out to be constant functions, not dependent on r: 

G(r) ~ G , (36) 

fi(r) ~ fj, Q , (37) 

w(r) =s wo- (38) 

The timelike component (/> t of the vector field has a negative exponential 
relationship to the radial distance: 



!5- 



— , (39) 

where Qs ac ts as a fifth- force charge. Note that in the large r limit, this result 
can also be obtained directly from Eq. (|19jl , after setting A' = B' = u>' — V$ (<j>) — 
Jt — and dropping the nonphysical positive exponential solution. 
The metric parameters A and B are in very close agreement with a Reissner- 
Nordstrom type solution for the spherically symmetric field of a charged mass: 



2G N M ojGqQ 



B(r)~l-=^^ + ^P±, (40) 



r 



A(r)~B(r)~\ (41) 

where Gn is Newton's gravitational constant, Q5 is the "fifth force" charge which 
is proportional to the source mass M. The emergence of the Reissner-Nordstrom 
solution, despite the fact that unlike electromagnetism, our vector field is massive, 
can be explained by the fact that at small values of r, mass contributions in the 
form of \ir become negligible, whereas at large r, Q\/r 2 rapidly approaches zero, 
and we are left with the Schwarzschild solution. 

In the next section, we will establish a relationship between Q5 and the source 
mass M. Using the values thus obtained, we find that the Reissner-Nordstrom 
solution is, in fact, a degenerate solution in which B{r) never vanishes, and which 
would normally include a naked singularity. However, at extremely low values 
of r, our numerical solution deviates from the Reissner-Nordstrom solution (see 
Figure [1]) 

The numerical solutions discussed here require that initial values be established 
for the functions involved. Specifically, we require, at some r = r, the values of G(f), 
G'(r), fi(f), n'(r), u{r), u'(f), A(r), B(f), B'(f), <f> t (r), and <j>' t {r). Note that A'{r) 
is absent from this list; the second derivative of A(r) does not appear in our set of 
equations. 

The values of A, B, and B' are fixed by the requirement that at large distance 
from a source, we must be able to mimic the Schwarzschild solution (albeit with 
a modified gravitational constant), and that at spatial infinity, the metric must be 
asymptotically Minkowskian. The vector field <fi must also vanish at infinity, which 
provides another boundary condition. Next, we assume that the values of G, /x, and 
ui are dependent on the source mass only, i.e., G' = fj/ = u>' = 0. In total, we have 
established seven initial values or boundary conditions through these assumptions. 

We seek the remaining initial conditions in the form of the fifth force charge Q$, 
and initial values of G = Go, fi = fJ-o, and lo — loq. We note that the basic properties 
of the numerical solution and the solution's stability are not affected by the values 
chosen for these parameters. However, their values must be chosen such that they 
correctly reflect specific physical situations. To determine these values, we now turn 
to the case of the point test particle. 
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5. Test particle equation of motion 

We begin by denning a test particle via its Lagrangian: 

£ T p = -m + aajq^^u 11 , (42) 

where m is the test particle mass, a is a factor representing the nonlinearity of the 
theory (to be determined later), ui is present as it determines the interaction strength, 
<75 is the test particle's fifth-force charge, and — dx^/ds is its four- velocity. 
We assume that the test particle charge is proportional to its mass: 

q 5 = Km, (43) 

with k constant and independent of m. This assumption implies that the fifth force 
charge q$ is not conserved, as mass is not conserved. This is the case in Maxwell-Proca 
theory, as V^J M ^ 0. 

From (|42p , the equation of motion is obtained p] [10] : 

m + r»pU a uA = -OLKumB» u u u . (44) 

That m cancels out of this equation is nothing less than a manifestation of the 
equivalence principle. 

Our numerical analysis of the spherically symmetric, static solution justifies 
the use of the Schwarzschild metric, provided that r ^> AiruiGoQ 2 . Further, for 
nonrelativistic motion, we can make the slow motion approximation ds ~ dt, yielding 
the equation 

•• J N GM 

r - + —3- = anuj(j)' t , (45) 

where r = dr/dt and Jn is the Newtonian angular momentum per unit mass. The 
second term on the left-hand side can be dropped at large r. Replacing cf>' t with the 
Yukawa solution (|39|) . and using ([43]) such that Q5 = kM, we get 

CM 

T 1-^(1 + F)e- F ■ (46) 



G 



Compatibility with Newton's equation of motion requires that when r <C /i 



-1 



from which 



r ~ -G N M/r 2 . Therefore, 

G (l-^)=G N , (47) 

k 2 = <LZ<™. (48) 
auj 

In particular, if we choose a such that 

G=(l + a)G N , (49) 

we get 



Given Q5 = kM, this allows us to write (f40|) in the form 
„, \ 1 2GnM G GnM 2 

B ( r > = 1 : — + — -2 — ■ ( 51 ) 
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Further, we also get 

Q 5 kuj = G N M. (52) 

Having obtained these results, we can now use the test particle to probe the spherically 
symmetric field of a point-like source. First, we note that for large r, A ~ B ~ 1, and 
the equations (|33|34|35|) for the scalar fields fi, oj, and G can be written as 

M" + V - ^ - %' ~ ^Gtfu? = 0, (53) 
r fi G 47T 

✓ ^-I^Otf + SgU ,54) 
G" + ?G'-^ + i(^-^G = 0. (55) 



Given fj,' = u>' = 0, ([55)) admits a constant solution. Conversely, holding 
G" = <p' t = 0, the equations for jj, and cj read 

r /i 47T 
2 , Gtffi 

oj H — cj H = 0. (57) 

r 07T 

If the dimensionless quantity G0j is small (which is the case if 4>t is calculated 
using (|39p for sufficiently large r), these equations uncouple from one another and 
admit approximately constant solutions, consistent with our numerical findings. 

We now choose to perturb these equations by introducing a constant radial 
component to the </> field, such that — ► — 4% < 0. For shorthand, we write 

4> 2 = 4>l - ti- (58) 
Still keeping the equations uncoupled, they are now solved by the following expressions: 



/ 47T 

/i = V^GW' (59) 

u = Ti-^Vr 2 , (60) 

48"7T 

2 

G = G °° (r + d)*- (61) 

How can these solutions be consistent with one another, and with our numerical 
results? First, we note that (1551) and (R)01) define a constant value for oj: 



Fixing oj in (p0|) also fixes r. We denote this value by f and note that it can be 
expressed as a function of M: 

f = KVM, (63) 

with K being a function of our chosen value of <f>: 



Fundamental parameter-free solutions in Modified Gravity 



10 



We choose <j) r to ensure that K remains constant, independent of M . Putting this 
into (f5"9")) we obtain 

^ = 7M' (65) 

with D given by 



Given [i oc M -1 / 2 , the value of p?Q§ is constant, since Q5 oc M. Given constant K 
in this means that the dimensionless quantity G(j> 2 must also remain constant; 

this, in turn, means a constant D. 

Conversely, from (frJTj) and ((63)) , we get 

G=G = Goo / ^ — , (67) 



with i? given by 



E 



Ci 
K 



(68) 



These two constants D and E are universal; their values are independent of M. 
Having established their constancy, we can choose to determine D and E directly 
from observation (This, , in turn, would fix the values of f and <f> r> although it is not, 
in fact, necessary to calculate explicitly the values of these quantities.) 

This result was obtained using a spherically symmetric, flat metric. Near a source, 
the flat metric is no longer a valid approximation; redoing the calculation using the 
Schwarzschild metric, we obtain 

G = G = G* + (Goo - G*) - j3- — , (69) 



(VM + Ef 

where agreement with Newtonian gravity demands G* = Gjv- 

These equations represent our results for the initial values of G, ui, and /x in our 
numerical solution: the numerical values of D and E, in turn, along with Goo, can be 
determined from observation. 

Our acceleration law therefore, can be written as 
G N M 



r 

with a, using (|49|) and (|69|) , given by 



[1 + a - a(l + fj,r)e-» r ] 



Ef 



Goo 

Gn 



1 



(70) 



(71) 



[VM 

and /i = no given by ([65]) . 

The acceleration law (170)) can also be recast in the commonly used Yukawa form 

GyM 



1 + ay 



,-r/A 



with the Yukawa parameters ay and A given by 

Gn 



G\ 



ay = 



1 + a Y ' 



(Goo - G N )M 



(G c 
= 1/ M = 



Gjv)M + Gjv(VM- 
M 



2' 



(72) 

(73) 
(74) 

(75) 
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We can also express the acceleration law (1701) as 

„ _ G eff M 

where the effective gravitational constant G c g is defined as 
G cff = G N [1 + a - a(l + He~' tr ] • 
The metric parameter B(r) can finally be written as 

B( r ) = 1 - 2GnM + ( 1 + a ) G NM 2 
This is the solution that was shown in Figure [T] 



(76) 
(77) 

(78) 



6. Cosmology 



In the case of a homogeneous, isotropic cosmology, using the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) line element, 



ds 2 = dt 2 - a 2 {t)[{\ - kr^dr 2 



l dVL 2 



the field equations assume the following form 

•■ 2 G 1 



jl + 3Hfi — 



/£ 

/'- 

G 



d> + 3Huj — —lu ■ 

Gr 



8?r 



G^ 



i 

47T 



gv + = o, 



(79) 

(80) 
(81) 
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G //i 



G + 3HG- -— + - ^ 
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8tt 



2 G 
3G (a 
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g Fg ~ y ^ 



G 



(82) 
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8ttGp 4tt / G 2 
3 T G 2 



+r GV, + 3 



8?r /Vb 



G 2 



A 



i? 
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4tt 
r G 



Gcu/j 2 ^ 



G' 



(83) 



4ttG 



(P + 3p) 



in & 
3 I G 2 



2 J. d V <t> A T 

tJH 0o - uj—— = At: J , 
o<Po 



A G G G 2 
+ 3 +H 2G + 2G ~ G 2 ' 



= J*. 



(84) 
(85) 



where H — a/ a is the Hubble expansion rate. 

Just as we were able to do so in the spherically symmetric case, it is possible to 
obtain an exact numerical solution to this set of equations using numerical methods. 
While noting that the stability or qualitative properties of the solution are not sensitive 
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Figure 2. The MOG "bouncing" cosmology. The horizontal axis represents time, 
measured in Hubble units of H^ 1 . The solid (black) line is a/ao, the scale factor 
normalized to the present epoch. The dashed (red) line is G/Go- The inset shows 
details of the bounce, demonstrating that a smooth bounce occurs even as the 
matter density of the universe is more than 10 14 times its present value. 




/ 



-1.0- 



Figure 3. The effective equation of state w c g (dotted green line) and the 
deceleration parameter q (dashed blue line) in MOG cosmology. Inset shows 
q in the vicinity of the bounce. The behavior of w is similar, in accordance with 
q = (l+3w)/2. Also shown with a solid red line is the effective density normalized 
to the cosmic scale parameter, computed as (Go£t(t) 3 Peft)/(GjvioPo)- 



to specific parameter values, for consistency with cosmological physics, we choose the 



following initial conditions: 

t = 13.7 x 10 9 years, (86) 

a Q = ct , (87) 

Go = 6£rjv, (88) 

Mo =H , (89) 

= 1/V12, (90) 

0o - 0, (91) 

ao = Hqoq = 72ao km/Mpc s, (92) 

G = Ao = c0 - 0, (93) 
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V a = 0.07659537Gg/ig, (94) 

V„ =V u = V$ = 0, (95) 

A =0, (96) 

k = 0, (97) 

The values of Hq and the Hubble time to, as well as that of a, are the same as in 
standard cosmology. We used for luq the result from the previous section. For /zo, we 
simply assumed that it is the inverse of the Hubble scale a , which we believe to be a 
reasonable choice. 

More notable is our choice of Go. The value of 6Gn is motivated by two 
considerations. First, we note that given a baryonic matter content of Qb = Pb/ Pent — 
0.05 (where pb is the baryonic matter density, and p cr it = 3Hq/8ttG is the critical 
density required for a flat universe with k = 0), a gravitational constant six times 
that of Newton decreases p C rit by the same factor without altering pb and, therefore, 
boosts VLb to the value of ~ 0.3, consistent with dark matter observations. Second, we 
note that an effective gravitational constant of G e ff — 6Gjv at the Yukawa distance 
r = fx , after solving for a in l|77[). yields a ~ 19 and an effective gravitational 
constant of G c s — 20Gn at infinity. On superhorizon scales our solution is consistent 
with an Einstein-de Sitter cosmology with no dark matter or dark energy. (This would 
imply a vanishing Vq on superhorizon scales.) 

To carry out the solution, we assume a pressureless matter equation of state 
w = p/p = 0. We find that this results in a "bouncing" cosmology. The bounce can 
be fine-tuned by choosing an appropriate value for Vg, which we have done. This 
ensures that the universe reaches sufficient density in order to form a surface of last 
scattering. At this point, the equation of state w = is obviously no longer valid. 
However, introducing a mixed equation of state that also incorporates radiation does 
not alter substantially the qualitative features of our model. 

We emphasize that in our model, only baryonic matter is present, with a matter 
density of ~5% of the critical density. Nevertheless, the cosmology is flat, due in part 
to the increased value of the gravitational constant, and in part to the presence of the 
nonzero energy density associated with Vg- 

The values of fi and lu remain constant, while G has a different value prior to the 
bounce, undergoing a rapid (but smooth) change at the time of the bounce. 

We can define an effective equation of state by rewriting ([83]) and (|84|) as follows: 



R2+ k_ , = 8*gA* 
a A 3 

Ob 47r 

= - —(l + 3w cS )Gp cS , (99) 
a 3 

and solving for w e f{. We can also define the deceleration parameter as 

1=~. (100) 

We find that w e ff is confined between — 1 < w e s < 0; it reaches its maximum at 
the time of the bounce, going asymptotically to its minimum value at times both after 
and prior to the bounce. The deceleration parameter q is negative in the distant past 
(indicating slowing contraction) . It becomes minus infinity at the time of the bounce 
(corresponding to a = 0), and becomes negative again at the present epoch, indicating 
accelerating expansion (see Figure [3|) . 



Fundamental parameter-free solutions in Modified Gravity 

NGC1560 NGC2403 



14 




12 13 14 IS 16 17 18 19 20 

NGC1399 



5 10 15 20 25 30 35 



10 20 30 40 50 60 70 80 90 100 



Figure 4. Galaxy rotation curves for a small set of galaxies of varying size. 
Data points are marked as black dots, current rotational velocity estimates are 
represented by a solid (blue) curve, while the dashed (green) curve shows velocity 
estimates in accordance with our earlier work 1131 114| . Mass estimates are as in 
[13] [14], except for NGC 1399, for which a mass estimate of M = 5 X 10 11 Mq 
was used. Dotted (red) curve is the Newtonian rotational velocity estimate for 
these galaxies using the same mass estimates. Radial distances are measured in 
kpc, masses in Mq. 



That MOG admits a bouncing cosmology without exotic matter or quantum 
effects may appear surprising at first, but a form of this bouncing solution is already 
known. As the /x and u> fields remain constant in this solution, they can be eliminated. 
The remaining field equations are those of Brans-Dicke theory [TT] with wbd = — 87r. 
Past investigations (see, e.g., [T2]) show that, indeed, a bouncing cosmology in Brans- 
Dicke theory can be obtained when wbd < — 6. However, standard Brans-Dicke theory 
lacks a Vq potential, which is what gave us the ability to fine-tune the solution and 
achieve a high density universe at the time the bounce occurs. 

Brans-Dicke theory is known to be in conflict with solar system observations, 
except when ubd is very large. However, we note that even though our cosmology is 
similar to Brans-Dicke cosmology, our spherically symmetric solution is not a Brans- 
Dicke solution. Here, the dominant contribution arises as a result of the M field, 
which is absent from Brans-Dicke theory. 
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7. Observations 

Unless one assumes that a massive dark matter halo is present, a typical spiral galaxy 
is dominated in mass by the central bulge. The motion of stars in the outer reaches 
of a galaxy can, therefore, be well approximated by the equations of motion in a 
static, spherically symmetric vacuum field. Indeed, our experience shows that the 
flat rotation curves of galaxies provide a sensitive test to determine the values of the 
constants D and E. In particular, it is easy to see that our results so far are compatible 
with the Tully-Fisher law [T5] . 

Kepler's laws of orbital motion yield a relationship between circular orbital 
velocity v c at radius r from a mass M in the form 

v 2 c GM , . 

7 = — ■ ( 101 ) 

In contrast, Tully and Fisher determined that for most galaxies, assuming that 
the brightness of a galaxy and its mass are correlated, the flat part of the rotation 
curve obeys the empirical relationship 

v n c tx M, (102) 

where 3 < n < 4. 

In our case, from (|6"3"]) and (|101[) we obtain for r ~ f: 

v\ cx -= = VM, (103) 
y/M 

corresponding to n — 4 in the Tully-Fisher relationship. 

Taking the next step, we have selected a small sample of galaxies that were studied 
earlier [13l [14] . An approximate fit to these galaxies yields the values 

D ~ 6250 M^kpc" 1 , (104) 

E ~ 25000 Mq /2 . (105) 

The galaxy rotation curves we obtain for galaxies of varying mass (including 
the recently studied galaxy NGC 1399 [16]) are in good agreement with these values 
(Figure [3} without dark matter. 

The galaxy rotation curves in Figure [4] were obtained modeling the galaxies as 
point masses, without benefiting from a core model, or the use of photometric data, 
as in the more extensive fit to galaxy rotation velocities [6]. Nevertheless, this exercise 
demonstrates that our newly established relationships between M, a, and \i not only 
satisfy the Tully-Fisher relationship in principle, but also offer good agreement with 
actual observations. 

Is it a mere coincidence that we have obtained good agreement with the Tully- 
Fisher relationship, or do our solutions work on other scales? Using the already 
established values of D and E, we attempted to reproduce the galaxy cluster mass 
profiles reported in 4 . Here, we run into a serious limitation: the solution we obtained 
is a spherically symmetric, static vacuum solution, and its applicability to an extended 
distribution of matter, such as that which one finds in a galaxy cluster, is limited. In 
a theory that offers linear behavior in the weak field limit, this is not necessarily a 
problem, as an extended distribution of matter can be well approximated by a large 
number of point masses. This is not so in our case: the gravitational field produced 
by two point sources is not simply the sum of their respective spherically symmetric 
static vacuum solutions. 
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Figure 5. A small sample of galaxy clusters studied in A ]. Thin (black) solid line 
is the mass profile estimate from |17| . Thick (blue) solid line is the mass profile 
estimated using our new results. Dashed (green) line is the result published in 
[1], while the dotted (red) line is the Newtonian mass profile estimate. Radial 
distances are measured in kpc, masses in Mq . 



This difficulty can only be addressed properly once appropriate interior 
solutions of our theory become known. Until then, our discussion is necessarily 
phcnomenological: the best we can achieve is a parameterization of the difference 
between any result derived from the spherically symmetric, static vacuum solution 
and a physically relevant interior solution. 

In U, the spherically symmetric, static vacuum solution was used successfully 
to model galaxy clusters. We are able to produce a comparable result while keeping 
the parameters D and E constant, by introducing an additional assumption: that the 
values of the MOG parameters and fi at some distance r from the center of a 
spherically symmetric distribution of matter are determined not only by the amount 
of matter contained within radius r, but by the amount of matter within radius r* . 
Figure \5\ shows the case of r* = 3r. 

As we have shown [3], our theory also appears consistent with key cosmological 
observations, including the luminosity-distance relationship of type la supernovae, the 
angular CMB power spectrum, and the galaxy-galaxy mass power spectrum (Figure[6|) . 

If the acceleration law (|16]) along with ([BIj]) and (|7ip is valid using the constants 
(I104|) and (|105[) . it should offer agreement with other observations. In particular, the 
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Figure 6. Cosmological observations and Modified Gravity (from [3]). Top 
panel: the luminosity-distance relationship of type la supernovae, with the MOG 
prediction shown with a thick (blue) line. (Thin (black) line is the ACDM 
prediction; dashed (red) line is a flat Einstein-de Sitter universe, while the 
horizontal access corresponds with an empty universe with no deceleration.) 
Bottom left: the angular CMB power spectrum showing good agreement between 
the MOG prediction and WMAP-3 and Boomerang data. Bottom right: After 
applying the appropriate window function, MOG (thick red line) shows agreement 
with the SDSS luminous red galaxy survey mass power spectrum, perhaps even 
superior to the ACDM prediction (thin blue line). 

following expression must be constant: 

[)?M = D 2 . (106) 

To verify the validity of this relationship, we have plotted M vs. ro = /i -1 in 
Figure [3 For the purposes of this plot, we used previously published results, while 
noting that our new calculations place dwarf galaxies, galaxies, and galaxy clusters by 
definition exactly on the line representing our prediction. This plot demonstrates the 
validity of Eq. (|106[) from the scales of star clusters to cosmological scales. 

The theory must also be consistent with observations or experiments performed 
within the solar system or in Earthbound laboratories. Several studies (see, e.g., |18j ) 
have placed stringent limits on Yukawa-like modifications of gravity based on planetary 
observations, radar and laser ranging, and other gravity experiments. However, our 
prediction of the absolute value of the ay parameter is very small when Ay is small. 
The latter is estimated at Ay ~ 0.16 pc (~ 5x 10 15 m, or about 33,000 AU) for the Sun, 
and Ay ~ 2.8 x 10~ 4 pc (~ 8.7 x 10 12 m, or ~ 58 AU) for the Earth. The corresponding 
values of \ay\ are \ay\ ~3x 10~ 8 and \ay\ ~ 9 x 10~ 14 , respectively, clearly not in 
contradiction with even the most accurate experiments to date (Figure [8]). 

While these observations do not prove that MOG is correct, they do suggest 
that MOG might be falsifiable. On the scale of Earth-based laboratories and in the 
solar system, experiments with ever greater precision might eventually rule out MOG. 
Beyond the solar system, as larger galaxy samples become available, the presence or 
absence of baryonic oscillations in the matter power spectrum may unambiguously 
decide in favor of modified gravity theories or dark matter. Confirmed detection 
of dark matter particles in deep space or in the laboratory would also be a strong 
indication against modified gravity. 
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Figure 7. The relationship fi 2 M = const, between mass M and the Yukawa- 
parameter rn = fi~ 1 across many orders of magnitude remains valid. The solid 
red line represents our theoretical prediction in accordance with H106I I. We are 
using cosmological data from [3]; galaxy cluster fits from [6]; galaxy and dwarf 
galaxy fits from 13 14] : satellite galaxy fits from :T ; and globular cluster fits from 
[9] . Note that the dwarf galaxy, galaxy, and galaxy cluster outliers are removed 
when these objects are recalculated using the results presented in this paper. 
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Figure 8. Predictions of the Yukawa-parameters from the MOG field equations 
are not in violation of solar system and laboratory constraints. Predicted values 
of A (horizontal axis, in m) vs. \cty\ are indicated by the solid red line. Plot 
adapted from 1181 . 



8. Conclusions 

In this paper, we have demonstrated how results of our Modified Gravity theory 
can be derived directly from the action principle, without resorting to the use 
of fitted parameters. After we fix the values of some integration constants from 
observations, no free adjustable parameters remain, yet the theory remains consistent 
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with observational data in the two cases that we examined: the vacuum solution 
of a static, spherically symmetric gravitational field, and a cosmological solution. 
These solutions were explored using numerical methods, avoiding the necessity to drop 
terms or make other simplifying assumptions in order to obtain an analytic solution. 
Further, the constraints used to compute the solutions are consistent to the extent 
that they overlap with one another, as discussed in Section [6l We view these results 
as encouraging, albeit not conclusive. The fact that at the level of the calculations 
presented here, our theory is not obviously falsified is an indication that it may be a 
worthwhile effort to pursue MOG further, for instance by obtaining interior solutions 
to the MOG field equations, or by developing tools to perform JV-body simulations to 
achieve more convincing results. 
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